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I. INTRODUCTION 

The calculation of finite temperature correlation functions is a long-standing problem in the theory of integrable 
models. Apart from the obvious conceptual importance of the problem there are many direct applications of the 
results to experiments on quasi-lD materials like KCuF 3 |ffl] or CuBenz which are described by integrable models. 
Finite temperature dynamical correlation functions in the systems have been measured by Neutron scattering and 
NMR and it is highly desirable to develop a method to calculate them exactly. 

Important progress was made during the eighties when it was realized that in gapless models Conformal Field 
Theory can be used to obtain the low-temperature asymptotics of dynamical correlation functions (see e.g. Q). 
However, models with a spectral gap as well as higher temperatures in gapless models remain outside the scope of 
the conformal approach. 

In a remarkable further development it became possible to determine the behaviour of static correlators at finite 
temperatures through an ingenious mapping of integrable d-dimensional quantum theories to d+l-dimcnsional inte- 
grable classical statistical models Q. However, dynamical correlation functions cannot presently be calculated by this 
approach. 

For integrable models with free fermionic spectra powerful methods to calculate finite temperature dynamical 
correlation functions have been available for some time |||| . Very recently the method of || was successfully extended 
to cases corresponding to interacting fermions @,^]. In particular, in (^] a formula describing the exponential decay of 
correlations in the delta-function Bose gas at finite temperatures was presented. Said formula is implicit in the sense 
that it is written in terms of solutions of certain nonlinear integral equations. The purpose of the present work is to 
analyze these integral equations by both analytical and numerical methods and present explicit expressions for the 
correlation lengths describing the decay of correlations. 

The outline of the paper is as follows. In section 1 we review some relevant facts on the delta-function Bose gas. In 
section 2 we study the special case of impenetrable bosons, in which particularly simple expressions for the correlation 
lengths are obtained at low temperatures. In section 3 we consider the general case of interacting bosons and we 
conclude in section 4. 



II. REVIEW OF THE 5-FUNCTION BOSE GAS 

The (^-function Bose gas is one of the paradigms of exactly solvable strongly correlated many-body problems in 
one spatial dimension. It describes N bosons interacting via a repulsive <5-function potential of strength c > 0. The 
Hamiltonian is 




At the special value c = oo the model describes noninteracting hard-core bosons and essential simplifications occur 
in the exact solution. We refer to this case as "impenetrable bosons" . The second-quantized form of the model is 
known as the Quantum Nonlinear Schrodinger equation and the Hamiltonian is expressed in terms of the canonical 
Bose field ip(x) as 
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H= dx(d x ^(x)d x ^(x)+cip Jf (x)^(x)i()(x)i()(x)) . (2) 



The model is solvable by Bethe Ansatz || and in what follows we recall some important ingredients of the exact 
solution. 

• Ground State and Excitations 

In momentum space a Pauli principle holds |L0| and consequently the zero temperature ground state is given 
by a filled Fermi sea of negative energy pseudoparticles. At c = oo this is simply the usual free-fermion ground 
state. The physics of the model is most conveniently described in terms of a rapidity variable A, which is related 
to the momentum k by 

fc(A) = A + [* 9{\-iM)f%(jjL)dii, (3) 

J-q 

where 9(X) = iln ( ^-a ) • density PtW °f pseudoparticles in the ground state is described by the integral 
equation 

PtW ^ fK{\n)p^)dn = i-. (4) 
Here q is the rapidity corresponding to the Fermi momentum kp and 

km- * + < 5 > 

Excitations over the ground state can be either particles or holes. The particle energy as a function of the 
rapidity is given by 



- ^ T K(\, M )£°(m) dp = X 2 -h, (6) 



where h > is the chemical potential. We note that the excitation energy vanishes on the Fermi surface 
e°(±g) = 0. The dependence of the integration boundary q on the chemical potential h is determined from the 
condition 

e°(±q) = . (7) 

In the impenetrable case the excitation spectrum is given in terms of free-fermionic particle-hole excitations. 
The Fermi velocity is defined as usual to be 



8e(X) 

Vf ~ 



dk{\) 



_ e°'(g) _ e Q, (g) ^ 
x=q k'(q) 2np° t (q) 1 1 



where the prime denotes differentiation with respect to A. The derivative of momentum with respect to rapidity 
on the Fermi surface is called dressed charge and is related to the density of states on the Fermi surface 

Z = 27rp° t (q)- (9) 



• Asymptotics of Correlation Functions 

The asymptotic behaviour of correlation functions at very low temperatures can be determined from the exact 
finite-size spectrum by conformal field theory techniques P,pd| . The result is found to be 

\x-v F t\ \x + v F t\\. (10) 
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The conformal dimensions A* are related to the dressed charge by 
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2A+ = 2A" = 



(11) 



In the framework of the Luttinger liquid approach these results were obtained by F. D. M. Haldane in |Tq] . 

In Q the following formula describing the exponential decay of correlations at any temperature was derived 
from the determinant approach to quantum correlation functions [^) 



^(0,0)^( X) t)) T ^ expj i- 



dX 



2np t (X) 



\x - v(X)t\ In 



e T 



e T 



exp (x) 



(12) 



The functions e(A), pt(X) and w(A) are the finite-temperature equivalents of the dressed energy, pse udoparticle 
density and Fermi velocity defined above. They are solutions of the following integral equations |13| 



e(A) = A 2 -h - — I dfi K{X, fj) In ( 1 + 



2tt 



Pt(A) 

The velocity v(X) is given by 



1 1 

2tt 2tt 



Km Ptifj) 

1 + e t 



v(X) 



1 <9e(A) 



27rp t (A) dX ' 

Note that for T — ► these equations reduce to (||) , (||) and (||) respectively. 



(13) 



(14) 



III. IMPENETRABLE BOSONS 



We first investigate the simpler case of impenetrable bosons. In the limit c — > oo (JL2J) simplifies to 



^(0,0)^(x,i)) T ^ exp<j — / dX \x-2Xt\ In 



A -h 

e t — 1 



e t 



(15) 



and we are left with performing a single integral. Expanding the exponentially decaying part of the integrand as 



In 



e t — l 



e t 



„ 1 
-2 > exp 

^2n-l 



T 



(2n- 1) 



(16) 



we obtain an expansion of (|15| ) in terms of error functions and exponential functions, which in turn yield an asymptotic 
low-temperature series on x m powers of T. The Fermi velocity is vp = 2yh, and we distinguish two cases. 

• Space-like region x/t > vp: 
After some elementary calculations we find the following expansion for \ ( as defined in (|l2|)) 

2m 

-2-2m\ 



^^yWflV (i_2- 

ttvf £-< n (2m ! \4hJ K 



l ) C(2m + 2) , 



(17) 



where £(ir) is the Riemann zeta-function. Note that this expansion breaks down as we approach x/t — > vp. The 
first few terms are 



xTn 
2v F 



5) WSl + 5049 ^) +0(2*) 



(18) 
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The inverse correlation length thus has an asymptotic power series expansion in odd powers of temperature for 
T — > 0. We also see that in the space-like regime there is no exponential decay with respect to the time t at low 
temperatures as the ^dependence only enters in terms such as exp(const/T) which are smaller than any power 
of the temperature T . 

• Time- like region x/t < vf: 
In the time-like regime we find that there are no power-law corrections to the conformal result 

X = ~ ^ exp (-[h - x 2 /4* 2 ]/T) + . . . (19) 

Thus, up to exponentially small corrections the conformal result is exact in the time-like regime and there is 
exponential decay of correlations with respect to t only at small temperatures. 



For general finite temperatures one needs to resort to numerical integration of . If we define a correlation length 
((T,h,i) by 

*-^-'*&krr«'*- x) wxs- (20) 

where <d(x) is the Heaviside theta-function, we can study the dependence of £ on T and the "direction" t/x by 
determining x numerically. From the above low-temperature analysis we already know that at very low temperatures 
£ becomes independent of t/x. In Fig [l] we plot the correlation length £ as a function of temperature for various 
values of t/x for the special values of chemical potential h = 1. 




FIG. 1. Inverse correlation length as a function of temperature for impenetrable bosons at h — 1. 

We clearly see the deviations from linear-T behaviour as temperature increases. The dependence on t/x becomes 
more pronounced at higher temperatures. In the regions tvF x and x <^ vpt the dependence on t/x disappears 
as expected. We also see that the exponential decay of correlations is fastest in the direction of the "light-cone" 

x/t = Vp. 
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IV. INTERACTING BOSONS 



We now consider the case of finite coupling constant c < oo. In order to derive an analytical low-temperature 
expansion of the correlation length we first need to perform an (asymptotic) expansion of e(A), v(X) and p*(A) in 
powers of T. 



A. Low temperature expansion of integral equations 



The expansion of e(A) in powers of T is readily established along the lines of jF4|. We first note that the function 
e(A) has precisely one zero for A > jH). We denote the corresponding value of the spectral parameter by qr- As e(A) 
is a symmetric function this then implies that e(±qr) = 0. From (O) and (ph we find 



e{\) - e°(A) - -L J dp K(X, p) [ E Qj) - e°(/*)] 
~ ~2^ 



dp K(X, p) In 1 + exp 



e(p) 



T 



T 
2^ 



d,K(X,,)ln + gLr) , (21) 
it l+exp(-|-^|) 



The last two terms in (^l|) are easily evaluated by expanding the integrand around ztqx- The leading contribution to 
the first term in (|2l| ) also comes from the regions p = ±qr and we again expand the integrand around these points 
and perform the resulting integrals. This gives 



s(X)-s°(X) 



2vr 



-Q 
7T 



12e'(q) 



[K(X,q) + K(-X,q)]T 2 + 



e'(gr) 

47T 



[K{X, q T ) + K (-X, q T )] (q T - qf + 0(T 3 ) . 



(22) 



Setting A = q in ( pl| ) we find that q — qx — 0(T 2 ) so that the last term in (21) does not contribute to 0(T 2 ). Putting 
everything together we then have 



s(A)= e °(A) 



1 



u(A) = K(X, q) + K(-X, q) + — / d/x X(A, /x) 



It immediately follows that 



q T = q 



12 [ e o'( g )]2 



(23) 



(24) 



The finite-temperature corrections to the pseudoparticle density T 2 Ap(p) = pt(p) — p^p) can be determined in a 
similar way 



Ap(A) 



12 [eO'( g )]s 



e°"(g) , p?'(g) 



2n e°'(q) p° t (q) 



u(X) - w(X) } + 0{T) 



(25) 



with 



w(X) = K'(X, q) + K'(-X, q) + — dp K(X, p) w(p) 

Zir 



where K'(X, q) — cLftT (A, q)/dX. Finally, the low-temperature expansion for v{X) is found to be 



v(X) = 





2^p°(A) 24 £ °'(g) p?(A) * 2n[p t (X)? 



1 u'(X) t2 e°'(X)Ap(X) T2 + 0{Jl3) 



(26) 



(27) 
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B. Asymptotics of the correlator at small temperatures 



Having obtained the leading corrections (in temperature) of e(A), PtM and U (A) we are now in a position to 
determine the leading temperature correction to the conformal result (|10|) for general values of the coupling c. In 
order to do so we need to expand the integral in (J12|), which is of the form 



1 

2^ 



dX /(A) In 



e T 



e T 



(28) 



where /(A) = ^ ^^(a* • ^ nc l° a di n g contributions to /(T) come from the vicinity of the points ±qr, where e(A) 



(29) 



vanishes. Expanding /(A) in powers of T, using (J23|), ( J25| ) and (|27[), and integrating around ±<7t we find 

X = AT + BT 3 + 0(T 4 ) . 



We find that ^4 reproduces the "conformal" result (|l0|). There is no correction to order T 2 . B is a complicated 
expression as it reflects the operator content of the theory so that we give an explicit expression only in the case 
x vpt. 



B 



Tr 3 Z'e "(q) 



3 71' 



■K A Z 



8Z 2 (e"'(q)) 4 24Z 3 (e°'((7)) 3 48Z 2 (e°'(q)) 3 8Z(e°'(q)) 5 
TT 2 u(q)e°"(q) _ ir 2 u'(g) | n 2 u(q)Z' ir 3 e°'"{q) 7r 2 Ap(g) 
24Z(e°\q)) 4 24Z(s°'(q)) 3 + 2AZ 2 (e ' (q)) 3 + 24Ze°'(q) + Z 2 e°\q) ' 



(30) 



where Z' = 2tt p° t '(q) and Z" = 2np° t "(q). 

At higher temperatures we again resort to a numerical solution of the relevant integral equations and integrals. In 
Fig. [2] we plot the inverse correlation length (|2(]) as a function of temperature for c = 1 and h = 1 . 




FIG. 2. Inverse correlation length as a function of temperature at c = 1, h = 1. 



G 



We obtain a qualitatively similar picture to the impenetrable case. Once again the exponential decay of correlations 
is fastest in the direction x — Dpt. However, the dependence on t/x is more pronounced than in the impenetrable 
case. This particularly clear in the space-like regime. 
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FIG. 3. Inverse correlation length for several values of c as a function of temperature at h = 1. 



V. CONCLUSIONS 



We have studied the rate of exponential decay of finite-temperature dynamical correlation functions of local fields 
in the (^-function Bose gas. We have explicitly evaluated corrections to conformal behaviour at low temperatures. For 
impenetrable bosons the results are qualitatively the same, and in addition we are able to get all the terms in an 
asymptotic series for the low-temperature behaviour of the rate of decay. 
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